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Counting the number of all the matchings on a bipartite graph has been trans- 
formed into calculating the permanent of a matrix obtained from the extended bi- 
partite graph by Yan Huo, and Rasmussen presents a simple approach (RM) to 
approximate the permanent, which just yields a critical ratio 0(na;(n)) for almost 
all the 0-1 matrices, provided it's a simple promising practical way to compute this 
^P-complete problem. In this paper, the performance of this method will be shown 
when it's applied to compute all the matchings based on that transformation. The 
critical ratio will be proved to be very large with a certain probability, owning an 
increasing factor larger than any polynomial of n even in the sense for almost all 
the 0-1 matrices. Hence, RM fails to work well when counting all the matchings via 
computing the permanent of the matrix. In other words, we must carefully utilize 
the known methods of estimating the permanent to count all the matchings through 
that transformation. 
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I. INTRODUCTION 



Let G = (y, E) be a bipartite graph, where V = ViU V2 i?> the set of vertices and 
-E C Vi X V2 is the set of edges. In the following sections we suppose #Vi = = n if 
there's no special illustration. A set of edges S <Z E is called a matching if no two distinct 
edges 61,62 G S contain a common vertex. S is called a k- matching if = k. In special 
case, S is called a perfect matching if k = n. Let Sk be the set of k-matching in G and 
A{G) be the set of all the k-matching, k = 0,1, . . . ,n. For the convenience of discussion, let 

n 

#S'o = 1, then the number of all the matchings in G is i^A{G) = J2 H^Sk- 
The permanent of a 0-1 A = aij,l < i, j < n is defined as 

n 

PeT{A) = Y,X{^^Mi) (1) 

TT 1=1 

where the sum is over all the permutations vr of [n] = {1,...,?t,}. It's well known that 
the permanent of an adjacent matrix of bipartite graph equals the number of its perfect 
matching. Let AM(G) denote the number of all the matchings in G, and A be adjacent 
matrix of G. [sj] has proved that 



AM{G) = ^per """" (2) 




where Inxn is n x n unit matrix, l„xn denotes nxn matrix with all the elements 1. This 
means in order to count the number of all the matchings of a bipartite graph with 2n vertices 
we only need to compute the permanent of a 2n x 2n corresponding matrix transformed 
from adjacent matrix. The computation of permanent has a long history and was shown 
to be T^P-complete in I2I]. Thus, in the past 20 years or so, many random algorithms 
have been developed to approximate the permanent, which can been divided at least four 
categories (3j]: elementary recursive algorithms(the original one is Rasmussen method(RM)) 
4]; reductions to determinants [sl, Q, 9, 11 1; iterative balancing 13]; and Markov chain 



Monte Carlo 



13 



All these methods try to find a fully-polynomial randomized 
approximation scheme fpras for computing the permanent, fpras is such a scheme which, 
when given e and inputs matrix A, outputs a estimator(usually a unbiased estimator)F of 



the permanent such that 



3 

Pr((l - e)per{A) < F < (1 + £)per{A)) > - (3) 



and runs in polynomial time in n and e"^, here 3/4 may be boosted to 1 — 5 for any 
desired 5 > by running the algorithm 0{log{5^^)) and taking the median of the trials 
jlol |. Then a straightforward application of Chebychev's inequality shows that running the 
algorithm 0{^^^e~'^) times and taking the mean of the results can make the probability 
more than 3/4(e.g. running A^^^e'"^ times). Hence, if the critical ratio -f^^ is bounded 
by a polynomial of inputs A, we'll get an fpras for the permanent of A. Another modified 
scheme called fpras for almost all inputs means: choose a matrix from A{n, l/2){A{n, 1/2) 
denotes a probability space of n x n 0-1 matrices where each entry is chosen to be 1 or 
with the same probability 1/2), or equivalently choose a matrix u.a.r. from A{n) {A{n) 
represents the set of n x n 0-1 matrices), and the following 

Pr(critical ratio of A is bounded by a polynomial of the input A)= 1 — o(l) as n — > oo 
holds. (Note that this is a much weaker requirement than that of an fpras). If a proposition P 
relating to n satisfies Pr(P is true)= 1 — o(n), we say P holds whp(whp is the abbreviation 
of "with high probability"). Thus, that there is an fpras for almost all the matrix means the 
critical ratio of A is bounded by a polynomial of the input A whp. A exciting result, that 
Markov Chain approach led to the first fpras for the permanent of any 0-1 matrix(actually 
of any matrix with nonnegative entry) was shown by [3]. However, its high exponent of 
polynomial running time makes it difficult to be a practical method to approximate the 
permanent. RM and reductions to determinants seem to be two practical approaches esti- 
mating permanent due to their simply feasibility, and both of them have been proved to be 
an fpras for almost all the 0-1 matrices, besides, Sj promises a good prospect on computing 
permanent via Clifford algebra if some difficulties can be conquered. RM also has developed 
to be a kind of approaches called sequential importance sampling way, which is widely used 
in statistical physics, see[l^. 

In this paper, we'll, by RM, compute the number of all the matchings based on the 
above transformation and give its performance theoretically, say, an analysis of critical 
ratio in the sense "for almost all the 0-1 matrix" of that matrix with a special structure. In 
section II, A new alternative estimator operating directly on the adjacent matrix without any 
transformation will be presented and proved to be equivalent to approximation performing 
on the transformed matrix by RM. In section HI, a low bound of the critical ratio for almost 
all the matrices will be presented, which is larger than any polynomial of n with a certain 
probability. Hence, RM does not perform well in computing the number of all the matchings 



as in computing the number of perfect matching. In section IV we'll propose some analytic 
results w.r.t. the expectation and variance of the number of all the matchings of a matrix 
selected u.a.r from Q{m,n){Q{m,n) denotes the set of bipartite graph with ^Vi = ^^¥2 = n 
as its vertices and exact m edges). These results seem likely to contribute to the upper 
bound of critical ratio for almost all matrices, but the calculations are more arduous and 
will be left for latter paper. 

II. AN EQUIVALENT ESTIMATOR 

All the notations have the same meanings as those in the previous section without 
special illustration. Let A an n x n 0-1 matrix be an adjacent matrix of a bipartite graph 
G — (V, E), {V — Vi[J V2). Set Ya a random variable. Then RM can be stated as follows: 

inputs: A an n x n 0-1 matrix; 

outputs: Ya the estimator of permanent A; 

if n=0; then 

Ya = 1 
else 

W^{j: ai, = 1} 
if = then 
Ya^O 
else 

Choose J u.a.r. from W 
Ya^\W\Yu 

Yij denotes the submatrix obtained from A by removing the 1st row and the jth 
column. Note this heuristic idea comes from the Laplace's expansion. Our following 
algorithm(for easy discussion, call it AMM) is also inspired by another expansion, we first 

presents our algorithm for the number of all the matchings, and then give the explanation 
and proof of equivalence between AMM and RM on the transformed matrix: 



inputs: A an n x n 0-1 adjacent matrix of G; 



outputs: Ya the estimator of the number of all the matchings of G; 
if n=0; then 

Ya = 1 
else 

W^{j: a^j = 1} U{0} 
Choose J u.a.r. from W 

Ya = \W\Yu 

Yio denotes a submatrix of A by removing the 1st row(of course, it's not necessarily a square 
matrix). Define a new terminology AM on the matrix, let B = {bij, 1 < i < m, 1 < j < n} 
an m X n matrix, m < n. let AM{0) — 1, by induction on m. 

n 

AM{B) := AM(Sio) + ^ feijSi,- (4) 

i=i 

Then we have the following theorem. 

Theorem 1. Let A be an n x n adjacent matrix of a bipartite graph G, Then 
AM(A)is the number of all the matchings of G. 

Proof: It's easy to check, when k > 1, the number of k-matching of G equals 
S ■ ■O'ikMik)^ where ii < i2 • • • < ik chosen from {l,2,---,n}, tt de- 

ii,— ,ik 1" 

notes the permutation of{ii, ^2, • • • , ife}. Thus, the number of all the matchings 

n 

is^ ^ X] ^n,7r(ii) ■ • ■ ^^ifc,7r(ijt) + 1; where 1 denotes the number of 0-matching. 

A;=l jl<---<ifcC{l,--- ,n} TT 

Note that if the AM(A) is written in terms of sum of elements of the matrix A, then it's 

n 

clearly to see AM (A) = E E aiiMn) ' ' ' ^ikMik) + 1-'^ 

k=lh<-<ikC{l,- ,n} TT 

Corollary 1. Let A — {oyT < i,j < n} be an n x n 0-1 matrix and Ya is obtained 
by above AMM. Then Ya is unbiased for AM(A), E(Ya) = AM(A) 

Proof: We prove for any m x n 0-1 matrix A, 1 < i < m, 1 < j < n}, which will be widely 
used in the following proves. AMM is unbiased for AM (A). For any fixed n, by induction 
on m, k=0,Vl < / < n, V a /c x / 0-1 matrix A, the equation E{Ya) = AM{A) is trivial. 
Now suppose yk < m,k < I < n, a, k X I 0-1 matrix A has E{Ya) — AM (A). Then when 
k — m, let \W\ — q, we have 



E{YA) = Y,E{YA\J = j)Pr{J = j) 

jaw 

jew 
jew 
jew 



- AM (A). 

□ 



Another simple corollary can also be obtained. To estimate the number of all the 

/ 



A Inxn 



matching in G, by RM operating on B 

1 n X ri 1 n X n 

operating on A by AMM, in precise words, which can be stated as follows. 



divided by n\ is equivalent to 



Corollary 2. Let Xa be the output of RM operating on A , Yb be the output of 
AMM operating on transformed matrix B divided by n\. Then Xa and Yb has the same 
distribution. 

Proof: Note that by RM after n-th step operating on B, Yb — Sn* ^inxn/^^'j where -S'^ is a 
number obtained from the first n steps, and obviously Yi^^^ = n\. Hence, we have Yb = Sn- 
The same distribution of Sn and Xa can be verified step by step.D 



A L 



nxn 



CorollaryS. AM (A) = ^per _ 

Inxn 1 n X n 

Proof: This is a direct deduction of corollary2. Let Xa be the output of RM operating on 
A , Yb he the output of AMM operating on transformed matrix B divided by n\. 

( 



AM{A) = E{Xa) = ECYb) = ^per 



^ ^nxn 



^ Inxn ^nxn 



□ 

So in the following section, we'll use AMM to compute all the matchings instead of RM 
since some methodologies similar to Rasmussen can be utilized. Another small advantage 



by AMM is that the critical ratio is smaller than that directly obtained from RM. The 
critical ratio by RM would be (2n)!, see Theorem 2.2^], while the critical ratio by AMM 
would be (n + 1)". 

Theorem2. Let A = {aij,l < i,j < n} he a.n n x n adjacent matrix of a bipartite 
graph G, and let Xa be the output of AMM. Then < {n + 1)". Generally, Let A be 

an m X n 0-1 matrix, m <n. be the output of AMM. Then ^^'^x^^ < + l)™ 
Proof: Induction on m. For any fixed n. k = 0,V1 <l<n, "iakxl 0-1 matrix A, the 
inequation is trivial. In the case k = m, let \W\ = q, we have 

E{Xl) = J2E{Xl\J = j)Pr{J = j) 
jew 

= J2e{xI,)i 

jaw 

<Y,E{XM,)\n + ir-\ 

j&W 

<(Ei?(X^J)^(n + ir-ig 

j<m 

= E{XA)\n^\r 

□ 

III. A LOWER BOUND OF CRITICAL RATIO FOR ALMOST ALL THE 

MATRICES 

Rasmussen shows that although the critical ratio of RM is factorial in n, it does indeed 
provide an fpras for almost all the matrix. However, the similar result can not be anticipated 
when computing all the matchings by RM. In fact the critical ratio for almost all the matrix 
would be more than rV^I'^~^ with a certain probability. To prove this, we need to define 
some new denotations. Since there're two probability spaces, we use the subscript a denote 
the calculus w.r.t. the probability space the algorithm lies in, say, coin-tosses, and subscript 
A represent the calculus w.r.t. the space probability the random matrices lie in. ^(m,n,p) 
denotes the probability space of all m x n 0-1 random matrices where each entry is chosen 



to be 1 with probability p, and^(m, n) denotes the set of all m x n 0-1 matrices . 

To obtain the mean and variance of the output of AMM on average under probabihty 
measure Pr^, we need the following lemma. 



Lemmal Let f{m, n) defined as f{m, n) — anf{m — 1, n) +Cnf{m — 1, n — 1), where m <n 
are two nonnegative integers, a„ and c„ are two infinite positive series w.r.t. n. And V 
0<l<n, f{0,l) = 1. Then 

m 

f{m,n) = J2 E c„ ■ • • Cn-fc+i< ■ • • + < 

k=l so+siH \-Sk=m—k 

so,---Sk>0 

Proof: By induction on m. Obviously, the case p=0 is trivial. Suppose when p < m — 1 

p 

V p < / < n, f{p, I) — E Q • • • Q-ifc+iaf° • • • a^^j^ + of holds, then when p=m, 

k=l so+siH hsk=p-k 

we have 

m— 1 

a„/(m - 1, n) = ^ c„ • • • c„_jk+i<°+^ • • • a'^^^ + al^ 

k=l so+siH \-Sf^=m—\—k 

m—l 

k=l so+siH \-Sk=Tn—k 

so>l 

and 

m—l 



... , „m-l 

"n-l-fc "-n^n-l 



Cnf{m-l,n-l) = ^ ^ c„ • • • c„_fca^°_i • ■ ■ < 

fe=l so+siH \-Sh=m—l—k 

m—l 

= C„ ■ • • Cn-kO'n-l ■ ■ ■ a^J^-fe + c„a„_i 

fc=l hsfe+i=ni— 1— fc 

m 

= ^ ^ c„ • • • c„_fc+ia^^_i • • • a^^_j^. + c„a^_]J^ 

k=2 s\+S2-\ hSfc=m— fe 

m 

— c„ • • • c„_fc+ia^Li • • • a^n_^. 

k=l s\+S2-\ hsfc=m— fe 



n—k 



k=l so+siH l-Sfc=TO— fe 

so=0 



Prom the above two equation, there holds 



f{m, n) = a„/(m -l,n) + c„f{m - 1, n - 1) 



k=l so+siH \-Sk=m—k 

so,-Sk>0 

The previous n can be replaced by any I, where m < I < n H 

Using lemmal we can easily obtain two following Theorems. 

Theorems. Choose A^^^ u.a.r. from A{m,n), m < n, or equivalently let ^^.n 
from ^(m, n, 1/2). Then 



EA(AM(Am,n))^yC^ — 



pk 

k=0 



where C^, = „ and = ^4at 

m k'.(in—ky. n (n— «)! 

Proof: Induction on m. The case p=0, £'^(AM(A)) = 1 is trivial. Suppose V p < m — 1, 
P<l<n E^iAM{Ap,i)) = E C^^fi = E 

k=0 k=0 

when p=m, \/ m < I < n, we have 



Using lemmal, here ai = 1, andc; = ^ then 



n 

77/ 

= £;^(^M(A^_i,0) + -EAiAM{Am-i,i-i) 

I 



Ea{AM{A^^i)) = J2 Y1 ■ ■ ■ ^'-^^+1 + ^ 



fe=l so+siH hsfe=m— A; 

™ pk 

fe=l so+siH |-Sfc=m— fe 

pk 



2k 

k=l 



EM pk 
2k 1^ 



k=0 



□ 

Theorem4 Choose A^n^^ u.a.r. from A{m,n), m < n, and let Xa^^^ be the output by 
AMM. Then 

™ pk 



2k 

k=0 



and 

771/ T~)lc T~)lc 

Ea{EMI^J) = J2 E + + 2 - 1)^^ ■■■{n + 2- ky^ 

k=0 so+siH \-Sk='m—k 

Proof: The first equation is is trivial since Ecr{X\^J — AM{Ajn,i)- For the second one, 
we use induction on m. The case p=0 is obvious. Suppose \fAp^i where < p < m — 1, 
p < I < n the second equation holds. When p — m, noting the fact M = — 1 is a 
binomial variable with parameter I and 1/2 (recall W^/{0} is the set of column indices with 
a 1 in the first row), then 
I 

Ea{E,{XIJ) = Y,Ea{E.{XIJ\M = q)PTA{M = q) 

q=0 
I 

= V EaUq + 1) V E^iX^i, )\M = g)Pr^(M = q) 
q=o jew 
I 

= J2 EaUq + l)E,{Xl_J + q{q + l)i?.(Xl_^_,_ J)Pr^(M = q) 

q=0 

= (Ea(M) + 1)Ea(E,{XI_J) + {Ea{M') + Ea{M))Ea{E,{XI_^^^J) 
- ('-^)Ea{K{XI^_J) + {^^)Ea{K{XI_^,J) 
Using lemmal, here a/ = and q = ^^^.Then 

fc=l so+siH \-Sk=m—k 

777 -r-^ ^ T^Jv 

= E^^ E {i + 2y%i + 2-iy^...{i + 2-ky^ 

k=0 so+si-\ l-Sfc=m— fe 

«o,---Sfe>0 

□ 

Theorems Choose An,n u.a.r. from A{n,ri), and let be the output by AMM. Then 

whp h{n) < E^{E,{Xa^J) < nh{n), where = M!^^^-^^, k* = [-1 + V^^l- 



where [*J denotes the largest integer no more than *. 
Proof: 

^ pk 

k=0 
" pn—k 

^n 2n-k 

k=0 

yn ( r, — 



2n Z_.(^_^)!(^!)2 

and let bk = j^^z$(j^: then ^ = set ^ > 1 we have A; < -1 + VS^TTS, thus, 

bk* = max bk- Thus, obviously 

A;=0,--- ,n 

□ 

Theorem6 Choose An,n u.a.r. from A{n,ri), and let be the output by AMM. Then 

whp 

Proof: Numerical experiment shows the above result, however the theoretical analysis 
seems so hard than until now I haven't thought out the way to show the comparably tight 
for since the order of J2 (n + 2)^o(n + 2 - l)^i • • • (n + 2 - A;)^Ms 

S() + Sl^ hSfe=n-fc 

too difficult to gain a good lower bound. The following bound is easy to check and the best 
one among methods I thought out, 

F(F(X^ ^wf- N)'(^ + 3)! 2^(A: + 2)^ 
J^A[^.[^A^,J) - (A;!)2(A; + 3)!(n - A;)! 

However it still can't reach the goal. Therefore, the proof of this theorem will be left for the 
future. 

Even if TheoremG has been proved, unfortunately, the critical ratio for almost all the 
matrices can not obtained from this theorem since two random variables are not independent. 
In order to accomplish the ultimate result, we need to calculate the ^)). Using 

the induction similar to theorem4, we can obtain the recursion of (recall M 

is a binomial variable with parameter n and |). 

E^{El{Xl^J) = 2{E^{M') + 2E^{M') + E^{M))E^{E,{Xl_jE,{Xl_^^^J) 
^E^(M')+2E^(M)+l)E^(El(Xl_J)+(E^(M')+2E^(M^^ 



Comparing £'_4(£'^(X^^ ^)) with E^{E^(X'^^J) and computing their ratio have to be 
done. Our main aim of doing this is to find the matrices satisfying E^{X\^J < 
Ej^{E'^{X\^^))g{n)^ where g{n) is a polynomial of n. However, the ratio of ^2 rx^'"'''')) 
so large that it can't accomphsh our goal. Thus we deduce our requirement whp to with 
a certain probability p > 0, and in our results p — \ — e where e is no more than 0.02. To 
prove the theorem, we need the following lemma, which will be proved in section IV. 
Lemma2 Let B{m,n) denote the set of all n x n 0-1 matrices with exact m I's, m » n. 
Choose B u.a.r. from B{m,n). Then 

n /^m—k 

r)2 



E{AM{B)) = J2iC'nn\ ^„ 

fe=0 



k 



and 

E{AM'^{B)) 



1 + o{l),n —>■ oo 



E^{AM{B)) 

Theorem? Choose u.a.r. from A{n,n), and let be the output by AMM. Then 



mix A. J- 2"^ 
where c is a constant no more 10, ands < 0.02. 

Proof: Prom lemma2 we know if we set m— (1/2 + e)v? and q — -g^- When n goes to 
infinity, noting k < n <^ 'm,n'^, there holds 

m(m - 1) . . . (m - k) 



and 



k-l 



ln{q) — ^^[Zn(m — i) — ln{'n? — i)] 

i=0 

= ^M5) + EMl-;j-)-Ml-;^)] 



= kln{-^)-J2[---,+0{—,)] 

i=0 

, , ,m, k(k-l) A 1 , /c^ , 

- kln{—) - -) + O — 



n 



3 



Thus, noting that km ^ < 2nm — 0{n^m ^) 

2 m 

_ (l/2 + £)n^ fc^ 1 1 

-( ^2 j expi 2^(i/2 + £)n2 ^ ^ ((1/2 + sV)^''^ 

<e-i(l/2 + £)'= 

Let B selected u.a.r. from B(m, n) Since g2^^|^|| = 1 + o(l), as n — > oo 

then Pr{AM{B) < |E(AM(5))) ^ 0, as n ^ oo. So, if m > (1/2 + s)^^ and £ < 0.02, we 

have whp 

E,{Xl) > EI{Xb) 
= AM\B) 

> {lE{AM{B))r 



r n /~im—k 
fc=0 



fc=0 

^E^^^ E (n + 2r(n + 2-ir...(n + 2-^r'= 

fe=0 so+siH hsfc=n— fe 



Ej,{E,{Xa^J). 



„2 



i=(l/2+e)n2 
3~ 



Noting Pr(>le U S(m,n))= 

m>(l/2+£)n2 

thus Pr(E.(Xl J > i?^(K(XA„,J)) > '^^^^y ■ 
Using Markov's inequahty. 



Pr(K(^l,J > nE^(K(^A„.J)) < - - 



then whp E„{X\^J < nEj^{E„{XA^,^)). Finally, we have 



■miXA^J - nE^{E,{XA„jy - 2-' 



Apply theorem6 to the above formula, we have 



k 



,2 



IV. THE NUMBER OF ALL THE MATCHINGS ON RANDOM GRAPH. 

In this section, we consider the expectation and variance of the number of all the 
matchings on G selected u.a.r. from Q{m,n). We have the following theorem. 

Theorems Choose G u.a.r. from Q {171,11), where Q{m,n) denotes the set of bipar- 
tite graph with #Vi = = n as its vertices and exact m edges, m ^ n, and let AM(G) 
denotes the number of all the matchings in G. Then we have 

n 

E{AM{G)) = Y,iC'j'k\E{XMik)) 

k=0 

and 

n k min{i,n—k) i—p 

k=0 1=0 p=0 j=0 

n k—1 min(i,n—k) i—p 

+ Y.Y.^^nfk\ Yl C:-kCrPL^-,pY^lp[Fn-A^-p-JmXMik^^-,)) 
k=l i=0 p=0 j=0 

where E{XMik)) = C'^l/CZ^ and = f:(-l)-C;Pr; 

r=0 

Proof: we'll use the methodology in [S]; Let M{k) be a k-matching on Vi + V2, For G G 
Q{Tn, n), define the random variable Xm{G) to be 1 if M{k) is contained in G, and otherwise 
0. The expectation and second moment of AM{G) is as follows. 

n n 

E{AM{G)) = EiJ2 E ^^m) = E E E{XMik)) 

k=0 M(k) k=0 M(k) 

and 

n n n 

E{AM\G)) = EiiY E ^^m)") = E E E ^(^M(.)^Mw) 

k=0 M{k) k=0 i=0 M(k)M'(i) 

where VO < A; < n, M{k) and M [k) range over all {C^fkl k-matching's on Vi + V2. Note 
that 

/^m~k 

E{XM{k)) 



G 



m 



The first equation follows quickly. For the second, in order to compute E(XM{k)X'^ii^i)), 
we have to calculate the number of pairs of M(k) and M' (i) as a function of the overlap 
j = \M(k) M'(z)|. For any fixed k, suppose i < k^ we need to compute the number of the 
pairs of M(k) and M' (i)^ where i = 0, ■ ■ ■ , fc, and M' (i) ranges over all (C^)^i! i-matching's 
on Vi + V2. The problem can be equivalently stated as follows: There're n different letters 
and n different envelopes. Among these letters, there're exact k{0 < k < n) labeled letters, 
each of which has only one 'mother envelope' among envelopes. Different labeled letters 
have different mother envelopes. We call a j-fit if there're exact j labeled letters put into 
its own mother envelope. Now choose i{0 < i < /e) letters from these n letters, then put 
them into i envelopes, and each letter can only be put into one envelope. V possible j, 
how many circumstances of j-fit are there? We can solve this problem like this: Suppose 
there're p letters unlabeled and i — p labeled letters among the selected letters, obviously, 
< p < min{n — k,i), the number of ways of choosing letters is C^^j^C]^^ . If the labeled 
letters has been laid, then the number of the ways of putting p unlabeled letters is -P^_(j_p). 
For any j{0 < j < i —p), there're C/_p ways putting exact j labeled letters in its own mother 
envelope. The last one we need to deal with is how many ways to put i — p — j labeled letters 
into n — j envelopes which contain all these i —_p — j letters' mother envelopes, satisfying 

0-fit. By the principle of inclusion-exclusion see[l|, we can easily obtain the number of the 

p 

ways is Fn-j{i — P — j), where = ^ (— l)''CpP^Z^. Noting that p ranges over to 

r=0 

min{i, n — k), and j ranges over to z — p, for each k and i < k. Then 

min(i,n~k) i—p 

M'(i) p=o i=o 



where E{XMik)) = CXl/Ci:. and = E (-l)^C;P„^i;. 

Consider, 



n n n k n—1 n 

EE E ^(^m(.)Xm«) = (EE+E E) E ^(^m(.)Xm«) 

k=0 i=0 M{k),M'{i) k=0 i=0 k=0 i=k+l M{k),M'{i) 

n k n—1 n 

= (EE+EE) E EiXMik)X'^,^) 

k=Q i=0 k=Qi=k+l M{k),M'{i) 

= (EE+EE) E ^(^M(.)x;,(,)) 

k=0 i=0 1=1 k=0 M{k),M'{i) 

-(EE+EE) E ^(^M(.)x;,,)) 

fe=0 i=0 k=l i=0 M{k),M'{i) 
n k n k—1 

fc=0 i=0 k=l 1=0 M{k)M'{i) 

- (E(^^n')^^!E + E(^n1'^!E) E ^(^M«Xm») 
fe=0 1=0 fe=l j=0 M'(i) 

Replace ^ E{XM{k)X'f^^^^) by 

M'(i) 

min{i,n—k) i—p 

E C^-kCr^Pn-i+p E Cl_p[Fn-j{i -p- j)]E{XM{k+i-j)) , then the second equation is 

p=0 j=0 

achieved. 
□ 



Remark: To complete the proof of theorem?, we also need to know whether the ratio 

E'\AM^(G)) to 1 as n goes to infinity, adding the condition such as m'^n~^ — > oo as n 

— >■ oo. We guess such a result is right, however the calculus seems very difficult. And this 
result also contributes to the upper bound of critical ratio for almost all the matrices. 
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